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ABSTRACT 

We present a detailed investigation of the incidence of circumnuclear dust 
structure in a large, well-matched sample of early-type galaxies with and with- 
out Active Galactic Nuclei (AGN). All 34 early-type AGN hosts in our sample 
have circumnuclear dust, while dust is only observed in 26% (nine) of a pair- 
matched sample of 34 early-type, inactive galaxies. This result demonstrates a 
strong correlation between the presence of circumnuclear dust and accretion onto 
the central, supermassive black hole in elliptical and lenticular galaxies. This cor- 
relation is not present at later Hubble types, where a sample of 31 active and 
31 inactive galaxies all contain circumnuclear dust. These archival, Hubble Space 
Telescope observations reveal a wide range of mostly chaotic dust morphologies. 
Current estimates suggest the dust settling or destruction time is on order of 
10^ years and therefore the presence of dust in ~ 50% of early-type galaxies 
requires frequent replenishment and similarly frequent fueling of their central, 
supermassive black holes. The observed dust could be internally-produced (via 
stellar winds) or externally-accreted, although there are observational challenges 
for both of these scenarios. Our analysis also reveals that approximately a third 
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of the early-type galaxies without circumnuclear dust have nuclear stellar disks. 
These nuclear stellar disks may provide a preferred kinematic axis to externally- 
accreted material and this material may in turn form new stars in these disks. 
The observed incidence of nuclear stellar disks and circumnuclear dust suggests 
that episodic replenishment of nuclear stellar disks occurs and is approximately 
concurrent with the fueling of the central AGN. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: structure — 
galaxies: elliptical and lenticular - dust, extinction — ISM: structure 

1. Introduction 

It is now widely accepted that most, if not all, galaxy bulges host a supermassive black 
hole (SMBH) at their center (Ferrarese & Merrit 2000; Gebhardt et al. 2000; Tremaine et 
al. 2002). However, only a fraction of these galaxies show evidence of significant accretion 
onto these central black holes in the form of Active Galactic Nuclei (AGN). The absence of 
significant accretion in many galaxies may be due to the lack of fuel (most likely in the form 
of cold interstellar gas), the lack of a fueling mechanism to remove angular momentum and 
drive this matter toward the center, or both. The most direct probe of the presence of fuel is 
the direct observation of cold interstellar gas, such as in the form of dust clouds in the central 
regions of galaxies. Cold dust is also dissipative and therefore sensitive to perturbations that 
may drive fuel toward the center of the galaxy. Thus the study of the dust distribution and 
its dynamical properties may provide valuable information about the powering, transport 
mechanisms, and timescales of nuclear activity. Identification of the origin of AGN fuel and 
the nature of the triggering mechanism(s) are two of the main unsolved questions in AGN 
and black hole research (e.g.. Martini 2004). 

A large fraction of all galaxies, including those with early-type morphology, have long 
been known to possess some form of dust structure (e.g., Kormendy & Djorgovski 1989). The 
best spatial resolution of the centers of nearby galaxies can be obtained with the Hubble Space 
Telescope {HST), which reveals the presence of cold interstellar clouds within parsecs of 
many nearby galactic nuclei. One of the first large HST studies of early-type galaxies by van 
Dokkum & Franx (1995) inferred that 78% of a sample of 64 galaxies contain circumnuclear 
dust. These authors also noted that the frequency of dust detections is higher in radio- 
detected galaxies (72%) than in galaxies without radio emission (33%), which suggested a 
correlation with nuclear activity. More recent HST studies have also found dust is common 
in the centers of early- type galaxies and a correlation between the presence of dust and line- 
emitting gas and/or AGN (Ravindranath et al. 2001; Tran et al. 2001; Lauer et al. 2005). 
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For later Hubble types, Malkan et al. (1998) also showed that dust structures were present 
in most HST WFPC2 images of a sample of 256 Seyfert and starburst galaxies; only a small 
fraction of the Scyfcrt galaxies did not appear to contain dust. Multicolor V— and i7— band 
observations of additional Seyferts further showed that circumnuclear dust was always found 
in the circumnuclear region of the predominantly spiral galaxy hosts of nearby Seyferts 2s 
(Regan & Mulchaey 1999; Martini & Poggc 1999) and Seyfert Is (Pogge & Martini 2002). 

A more robust way to address the relation between the presence of circumnuclear gas 
and/or dust and nuclear activity is to look for differences between AGN host galaxies and 
a control sample of inactive counterparts. In a recent paper, Xilouris & Papadakis (2002) 
searched for dust structures in HST images of a sample of 23 Seyfert and 35 inactive galaxies 
by fitting and subtracting elhptical isophotes. They found that all active galaxies possess 
some form of nuclear dust structure and, for their subsample of early-type galaxies, the 
Seyferts exhibit more nuclear structure than the inactive galaxies, although the opposite 
seems to be true for the later types. In a more recent paper. Martini et al. (2003) investigate 
the incidence of dust in a well-matched sample of 28 active and 28 inactive galaxies and 
found that all of the AGN possess dust structure, while seven of the inactive galaxies did 
not. 

Studies of the centers of early-type galaxies have also revealed that many contain disky 
isophotes, and in some cases these isophotes are most likely due to embedded stellar disks 
(e.g., Capaccioli et al. 1987; Bender et al. 1988; Kormendy & Djorgovski 1989; Scorza & 
Bender 1995; Seifert & Scorza 1996). Observations with HST have identified many stellar 
disks within the central kiloparsec (hereafter nuclear stellar disks) that are distinct from the 
outer disks present in lenticular galaxies (van den Bosch et al. 1998; Scorza & van den Bosch 
1998; MoreUi et al. 2004; Krajnovic & Jaffe 2004). With a large sample of 67 early-type 
galaxies. Rest et al. (2001) found nuclear stellar disks in as many as 51% of their sample, 
although many were misaligned and likely nuclear stellar bars. Given that disks are easiest to 
detect when nearly edge on, the observed fraction is consistent with the presence of nuclear 
stellar disks in all early- type galaxies (Rix & White 1990). 

In the present paper we address the relation between circumnuclear dust structures and 
the nuclear activity in galaxies, with a particular emphasis on early-type galaxies. Unlike 
most previous work, we use a well-defined and homogeneous input sample - the Palomar 
sample of Ho et al. (1995) - and have carefully drawn a pair-matched sample of 34 active 
and 34 inactive early-type galaxies (ellipticals and lenticulars), matched by the properties 
of their host galaxies, in order to obtain a robust measurement of the frequency of dust 
structures in active and inactive early-type galaxies. We have also identified and analyzed a 
pair-matched sample of 31 active and 31 inactive late-type galaxies. This paper is organized 
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as follows: In Section 2 we present the sample selection and the pair-matching technique; in 
Section 3 we present data reduction and analysis procedures; in Sections 4 and 5 we present 
our results and discussion, respectively, and finally in Section 6 we present our conclusions. 

2. Sample 

A key issue in the comparison of active and inactive galaxies is the selection of a well- 
defined control sample known to be inactive. This is difficult because identification of a 
low-luminosity AGN often requires less sensitive observations than confirmation that a low- 
luminosity AGN is not present. In order to avoid this problem, we have drawn our sample 
from the large, uniform, and sensitive Palomar Survey of Ho et al. (1995). The Palomar 
Survey contains spectra of the closest 486 bright galaxies in a set region of the sky and 
all of the nuclear spectra are classified as: absorption-line nuclei, HII, LINERs, transition 
objects, and Seyfert galaxies (Ho et al. 1997a). The Palomar Survey is considered to be the 
most complete and homogeneous representation of the nearby universe. Our first selection 
of the active galaxy sample comprised all Palomar Seyfert and LINER galaxies with avail- 
able broad-band HST WFPC2 images in the optical spectral region, excluding transition 
objects and those with uncertain classification. The sample of inactive galaxies comprised 
all Palomar galaxies classified as absorption-line and HII nuclei with available broad-band 
WFPC2 images in the optical spectral region. 

We then carefully selected a well-matched control sample from this sample of inactive 
galaxies through use of a pair-matching algorithm. This approach is similar to that adopted 
by Martini et al. (2003), who argue that the best way to match control samples is to identify 
a control galaxy pair for each active galaxy with similar values of all properties that may 
affect the identification of dust structure (e.g., morphology, distance, luminosity, inclination). 
More traditional and commonly employed techniques, in contrast, just match the mean or 
median value of each property between the target and control samples, or the distribution of 
each property individually, and may be susceptible to different correlations between various 
properties in the two samples. 

In order to make a robust comparison between active and inactive galaxies we have 
identified an inactive pair for each active galaxy with the following criteria: maximum dif- 
ference between absolute magnitudes of 1 mag; maximum difference between morphological 
T types of 1; maximum difference between galaxy inclinations of 15°; maximum difference 
in plate scale of 50%. To accomplish this we developed an algorithm that finds the largest 
possible number of pairs between the active and inactive galaxies given the above criteria. 
This yielded 26 early-type and 31 late-type pairs. These 57 pairs comprises our main sample 
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and we shall hereafter refer to it as the "matched sample" . In order to increase the early- 
type sample we then reran the pairing algorithm with the remaining early-type galaxies after 
relaxing the matching constraints to: maximum difference between absolute magnitudes of 
3 mag; maximum difference between T types of 3; maximum difference between galaxy in- 
clinations of 35°; maximum difference in plate scale of 80%. This resulted in 8 additional 
early-type pairs for a total of 34 matched pairs of early-type galaxies. This sample, com- 
bined with our 31 late-type pairs, we shall refer to as the "extended sample." Histograms 
of the Hubble types, distances, inclinations and absolute B magnitudes for the matched and 
extended sample are presented in Figures 1 and 2 respectively. 

These figures show there are similar distributions for each property in the active and 
inactive galaxies. Although this is expected for the matched sample, it also remains true for 
the extended sample, and demonstrates that the additional 8 early- type galaxy pairs in the 
extended sample are reasonably well-suited for comparison of the properties of active and 
inactive early-type galaxies. Relevant information for the sample galaxies is presented in 
Table 1 for the active galaxies and Table 2 for the inactive galaxies. Columns 5-9 of each of 
these tables list the galaxy T type, distance, absolute magnitude, inclination, and activity 
type from Ho et al. (1997a). 



3. Data Reduction and Analysis 

These data comprise EST WFPC2 images obtained both with the planetary and wide 
field cameras. Relevant data on the images, notably the camera, filter, and exposure times, 
are listed in columns 2-4 of Tables 1 and 2. All of the images were initially reduced by 
the HST archive's on-the-fiy reprocessing pipeline. The only additional processing was 
cosmic-ray removal with the appropriate IRAF tasks. Specifically, when two or more images 
were available we used the CRREJ task in the STSDAS package, otherwise we used the 
COSMICRAY task in the NOAO package. 

To identify circumnuclear dust structure in these images we used the structure map 
technique proposed by Pogge & Martini (2002). Structure maps enhance structure as fine 
as the scale of the point spread function in an image and are well suited to the identification 
of narrow dust lanes and emission-line gas in nearby galaxies. The technique is based on 
the Richardson-Lucy (R-L) image restoration algorithm, which uses the PSF to identify 
structure to enhance as part of the deconvolution process. Mathematically, a structure map 
S is defined as: 

-r^ P' (1) 
I®P ^ ' 
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where / is the original image, P is the instrument PSF (constructed with the TinyTim 
software of Krist & Hook 2004), P* is the transpose of the PSF and ® is the convolution 
operator (see Pogge & Martini 2002). 

The structure map technique recovers similar information to the more traditional meth- 
ods of contrast enhancement, such as fitting and subtracting elliptical isophotes and construc- 
tion of color maps, but offers several advantages. For example, the fitting and subtraction 
of ellipses can lead to artificial features in the presence of strong brightness discontinuities 
or isophotal twists in the images. In the case of color maps, there are simply not as many 
galaxies with observations in two filters as in one. Possible mismatches between the PSFs in 
the two bands could also introduce artificial features in color maps. These problems are all 
avoided with the use of structure maps, although some spurious features can still appear. For 
example, overexposed pixels can produce dark rings that mimic absorption features, while 
dead pixels can appear as false bright spots in the images. Nevertheless, such cases were 
easily identified and did not affect our analysis. After constructing structure maps, we then 
visually inspected each galaxy to determine if dust structures were present and measured 
the spatial extent of any dust. We also noted the presence of emission features that seem to 
be nuclear stellar disks (see discussion below) and measured the projected radial extent of 
these features. 



4. Results 

Structure maps for all our galaxies are shown in Figures 3 to 5. Figure 3 present 
structure maps of each early- type (T < 0), AGN host next to its inactive, control galaxy 
for the matched sample. Each galaxy pair has the active galaxy on the left and the inactive 
galaxy on the right. These two figures illustrate the most striking result of this paper, namely 
that all of the early-type active galaxies possess circumnuclear dust while only 27% (7 of 26) 
of the early- type inactive galaxies possess circumnuclear dust. Figure 4 demonstrates that 
this result also holds when we include the extended sample, specifically all of the early-type 
active galaxies have circumnuclear dust and only 26% (9 of 34) of the inactive galaxies in the 
extended sample have circumnuclear dust. This result also holds if we exclude lenticulars 
and only consider ellipticals (T = —5). Eight of eight (10 of 10) active ellipticals in our 
matched (extended) sample have dust, while only one of eight (two of 15) inactive ellipticals 
in our matched (extended) sample have circumnuclear dust. 

While previous investigators have noted an excess of circumnuclear dust in early-type, 
active galaxies, compared to early-type, inactive galaxies, ours is the first study to com- 
pare the frequency of circumnuclear dust between well-matched samples. In addition, both 
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our active and inactive, control samples have sensitive and uniform spectroscopy from the 
Palomar survey. These data are important to confirm both the presence of activity and its 
absence in the control sample. Some of the earliest evidence that early-type, active galaxies 
may more commonly possess dust came from studies of radio galaxies (Sadler & Gerhard 
1985; Marston 1988; JafTe et al. 1994). van Dokkum & Franx (1995) presented the first study 
of the incidence of dust in active and inactive early- type galaxies with HST and found dust 
in 72% of the radio galaxies, while only in 33% of the galaxies without radio emission. Other 
investigators (Verdoes Kleijn et al. 1999; Ravindranath et al. 2001) have similarly found 
circumnuclear dust in most active galaxies, while rarely found dust in inactive, early-type 
galaxies. The Ravindranath et al. (2001) sample is the most comparable to our own because 
that sample was also selected from the Palomar survey, although their sample is somewhat 
smaller (33 galaxies total), did not match active and inactive hosts, and not all galaxies 
possessed visible-wavelength HST images. We attribute our 100% detection rate of dust in 
early-type AGN hosts to the excellent sensitivity of the Palomar Survey to low-luminosity 
AGN. 

The late-type (T > 0) galaxies are shown in Figure 5. All of these galaxies display 
evidence for circumnuclear dust regardless of the presence or absence of detected nuclear 
activity. Previous studies of dust in galaxies with similar Hubble types had shown that 
circumnuclear dust is common in active, late- type galaxies (e.g., Malkan et al. 1998; Regan 
& Mulchaey 1999; Martini & Pogge 1999; Pogge & Martini 2002). Circumnuclear dust has 
also long been known to be quite common in late-type, inactive galaxies. With a similarly 
matched sample of 28 active and 28 inactive galaxies. Martini et al. (2003) found that all 
active galaxies had circumnuclear dust, although 25% (7 of 28) of the inactive galaxies 
did not show evidence for circumnuclear dust. Only three of the seven inactive galaxies 
in their sample have later Hubble type than SO, therefore their measured dust frequency is 
consistent with the presence of dust in all 31 of our late-type, inactive sample. The frequency 
of circumnuclear dust structure as a function of Hubble type for the matched and extended 
samples are shown in Figure 6. 

Most of the dust features in these early-type galaxies are dominated by dust lanes 
and in many cases the lanes trace some spiral structure that indicates the presence of a 
dusty disk. In total, at least half of the active early-type galaxies exhibit some evidence for 
dust disks and 25% of these have tightly wound spirals in regular disks. These later dust 
structures are interesting because tightly wound dust structures are generally observed in the 
most axisymmetric (i.e., unbarred) late- type galaxies and suggest the presence of long-lived 
structures (Peeples & Martini 2006). The presence of any dust structures in each galaxy 
is indicated in column 11 of Tables 1 and 2, where we report the projected radial extent 
of each dust feature in kiloparsecs. This radius was estimated by visual inspection of the 



-8- 



digital images. A x sign means that the radius of the structure is not clear, and may extend 
beyond the borders of the image. For the late-type galaxies we have not measured the radius 
of the dust structure as in most cases it does extend beyond the borders of the images. The 
dust disks and spirals have a mean size of 0.43±0.37kpc and 0.69±0.63kpc respectively 
for the early-type active galaxies. These structures therefore appear to be confined to the 
circumnuclear region or central kiloparsec, although these observations are less sensitive to 
dust on larger scales due to a decrease in signal-to- noise ratio. The position angles of the 
dusty disks in the early-types, if apparent, are generally aligned with the position angle 
of the large-scale isophotes. Some exceptions include orthogonal dust disks in NGC 2787, 
NGC 4111, and NGC 4143. 

In addition to circumnuclear dust, our structure maps also reveal the presence of nuclear 
stellar disks in many galaxies. We have found disks in at least 13 of the 34 (38%) early-type 
inactive galaxies, while we found nuclear stellar disks in only one of the 34 (3%) early-type 
active galaxies. Many of these nuclear stellar disks were identified in previous work, such as 
NGC 4570, NGC 4621 and NGC 5308 by Krajnovic & Jaffc (2004), which employed elliptical 
isophote fits to identify these disks. As we have identified a similar fraction of stellar disks, 
structure maps may be as sensitive to nuclear disks as isophote fitting. This is somewhat 
surprising, as structure maps arc best suited to the detection of features on the scale of 
the point spread function. The appearance of the stellar disks in the structure maps may 
therefore indicate that they are nearly edge on, as is the case for elliptical isophote fits. 
Consistent with this interpretation, nearly all of the nuclear stellar disks in our sample have 
very large axis ratios. The detectability of stellar disks in early-type galaxies was modeled 
in detail by Rix & White (1990), who found that the observed fraction is consistent with the 
presence of a nuclear stellar disk in every early-type galaxy. 

The apparent nuclear stellar disks in these galaxies could also be nuclear stellar bars 
(e.g., Erwin & Sparke 2002). To determine if this is the case we measured the position 
angle of each of the nuclear stellar disks and compared it with the position angle of the 
galaxy isophotes on large scales. In all cases the nuclear stellar disk was aligned with the 
position angle on larger scales, which suggests that they are not nuclear bars. Surprisingly, 
the position angles of the nuclear stellar disks in five out of five elliptical galaxies also agree 
with the position angle of the isophotes at large scales (NGC 821, NGC 3610, NGC 3377, 
NGC 4621, and NGC 7619). While some of these eUipticals are known to be disky (e.g., NGC 
3377, 4621), ellipticals are known to be triaxial and do not all rotate about their apparent 
minor axis (Ryden 1992). We return to this point in Section 5.2. 

In contrast to the early-type, inactive sample, we do not detect many nuclear stellar 
disks in the early-type active sample or in either of the late-type subsamples. The absence 
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of nuclear disks in late-type galaxies has already been reported in the work of Pizzella et 
al. (2002), who studied HST WFPC2 F606W images of a sample of 38 spiral galaxies 
and did not find nuclear disks in barred galaxies or galaxies of Hubble type later than Sb. 
These differences from our early-type, inactive sample arc most likely due to the substantial 
circumnuclear dust in the early-type active sample and the late-type galaxies. In almost 
all cases, the dust (and emission-line regions) present in these galaxies make it difficult to 
identify nuclear stellar disks, particularly if the dust and stellar disks have the same position 
angle. Some notable exceptions are NGC 4111, which has a prominent stellar disk almost 
perpendicular to a smaller dust disk, and NGC 4026, which has a nuclear stellar disk and 
weak dust features. If nuclear stellar disks are also present in the dusty, active early-types, 
the fraction may be higher in a study of near-infrared observations of early-type galaxies, 
although as Ravindranath et al. (2001) note the nuclear stellar disk fraction may also depend 
on the central surface brightness profile. 

A histogram of the presence of stellar disks as a function of Hubble type in our matched 
and extended samples is shown in Figure 7. The presence of nuclear stellar disks in our sample 
is identified in column 12 of Tables 1 and 2 and provides the radial extent of the structure 
in kiloparsecs. In Figure 8 we present a histogram showing the radial extent in kiloparsecs 
of the nuclear disks, which have a mean value of 0.24±0.31kpc and never extend beyond 
1.1 kpc. Therefore the bright stellar disks in the inactive galaxies are nuclear structures, 
similar or somewhat smaller in average spatial extent to the dust structures observed in the 
active galaxies. 



5. Discussion 

In the previous section we identified two significant differences between active and inac- 
tive early-type galaxies. First, all active early-type galaxies have circumnuclear dust, while 
dust is only present in 26% of inactive early-type galaxies. Second, 38% of the inactive 
early-type galaxies have nuclear stellar disks, while they are detectable in almost none of the 
active, early-type galaxies. In the subsections below we discuss these results in the context 
of AGN fueling, the origin of the dust and stellar disks, and present a simple evolutionary 
scenario for the centers of early-type galaxies. 
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5.1. Ubiquitous Dust in Early- Type AGN Hosts 

Dust is present in all of the active early-type galaxies, both in the matched and in the 
extended sample. In contrast, dust was only found in 27% (7/26) of the inactive, matched 
sample and 26% (9/34) of the extended, inactive sample. This strong correlation indicates 
that the presence of circumnuclear dust is a requirement for black hole accretion in early-type 
galaxies, while the presence of circumnuclear dust in some inactive galaxies indicates it is 
not a sufficient condition for activity. This strong correlation implies that the circumnuclear 
dust in these galaxies, most of which is hundreds of parsecs from the central black hole, is 
connected to or had a similar origin to the material currently fueling the active nucleus. 

As the circumnuclear dust observed in early-type, inactive galaxies almost always ex- 
tends to the unresolved nucleus (tens of parsecs), this suggests that when dust is present 
in early-type galaxies, at least some dust reaches the nucleus relatively quickly compared 
to its destruction time. The ratio of inactive to all early-type galaxies with circumnuclear 
dust in a representative (as opposed to matched) sample multiplied by the dust lifetime 
could provide an estimate of how long it takes circumnuclear fuel to reach the nucleus. If we 
take the early-type galaxies in the Palomar survey as representative of all early-types. Ho et 
al. (1997b) found that approximately 50% of all early-type galaxies are AGN. Our matched 
sample of early-type active and inactive galaxies is therefore likely representative of all early- 
type galaxies and we infer that the time for some fuel to accrete onto the central black hole 
is approximately 20% of the lifetime that dust is observable in the circumnuclear region. 
This percentage would be yet smaller if the dust in the circumnuclear region compromised 
the identification of some early-type galaxies as active. 

There are differing models for the origin and hfetime of dust in early-type galaxies. 
Unlike in late-type disks, there is not an obvious reservoir of dust on larger scales. The 
absence of such a dust reservoir indicates that the dust must be either created in situ or 
originate from outside of the galaxy. Stellar mass loss could produce the observed dust 
in early-type galaxies (Knapp et al. 1992). Athey et al. (2002) show that the mid-infrared 
dust emission from early-type galaxies is consistent with the expected stellar mass-loss rates, 
where the mass-loss rate is on order 1 MQyr~^ for a giant elliptical galaxy. The two difficulties 
with this process are that the dust should be produced in a uniform distribution that traces 
the stellar kinematics and the dust will gradually be destroyed via sputtering by the hot 
gas that pervades these galaxies. The dust distribution in a flattened disk is somewhat 
puzzling in the elliptical galaxies, such as NGC 4261, because they generally do not exhibit 
substantial rotation, although even initially slow rotation would lead to substantial rotation 
if the dust settled into the central kiloparsec from large scales. Mathews & Brighenti (2003) 
model the evolution of the dusty gas originating in stellar winds and found cooling dust-rich 
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gas could produce the observed circumnuclear dust, and in particular the observed clumping. 
However, even if dust cooling explains the observed dumpiness, stellar mass loss still can not 
explain why circumnuclear dust is only present in ~ 50% of all early-type galaxies. While 
dust destruction via collisions with hot gas will occur in on order 10® years (e.g., Mathews 
& Brighenti 2003), mass loss will continually supply new dust to the interstellar medium. 
This difficulty with internal dust production is also reflected in the absence of a correlation 
between far-infrared and visible-wavelength luminosities in early-type galaxies, which also 
suggests the dust may have an external origin (Temi et al. 2004). 

Kinematic observations of neutral and ionized gas in carly-typc galaxies have led many 
authors to conclude that this material has an external origin because the gas kinematics differ 
from the stellar kinematics (e.g. Bertola et al. 1984, 1992; Sarzi et al. 2006; Morganti et al. 
2006, and references therein). The circumnuclear morphology of the dust also provides some 
additional constraints on the origin of the dust. In particular, the chaotic appearance sug- 
gests an external origin because internally-created dust might be too uniformly distributed, 
particularly in lenticulars, if the cooling time discussed above is not shorter than the dynam- 
ical time. However, the circumnuclear dust observed in late-type galaxies exhibits a similar 
range in morphologies and that material is generally not held to have an external origin. 
The morphology of the dust can also provide some constraints on the lifetime of the dust. 
Lauer et al. (2005) discuss how the different dust morphologies can be viewed as a "setthng 
sequence," where dust that is not in a disk must settle toward the center of the galaxy in a 
dynamical time, or a few 10^ years (see also Tran et al. 2001; Verdoes Kleijn & de Zeeuw 
2005), while coherent, tightly wound dust disks must be at least a few dynamical times old. 
The age of such a dust disk would then be set by the dust destruction timescale mentioned 
above, or ~ 10® years. External dust input and a finite dust lifetime can plausibly explain 
why dust is not found in the centers of all early-type galaxies. 

More broadly, the relative number of galaxies with chaotic dust lanes, some nuclear dust 
spirals, and very coherent dust disks could provide constraints on the settling time of the 
dust if there is a natural progression toward more ordered dust as the dust settles toward 
the galaxy nucleus or is destroyed. If this model is correct, then the dust in inactive galaxies 
should be most chaotic or unsettled. Our data are consistent with this interpretation, with 
the exception of the tightly wound dust spiral in NGC 4526, although there are only a 
small number of tightly wound dust spirals in our full early-type galaxy sample. In fact, the 
relative scarcity of the tightly wound dust spirals presents a problem because the estimates 
of the dust settling time suggest that these structures should be longer-lived than the more 
commonly observed amorphous dust. As noted by Lauer et al. (2005), this suggests the 
total dust lifetime must be close to the dynamical time of several times lO'^ years and these 
structures can not be long lived. 



- 12 - 



Both the internal and external dust models rely on a dust destruction timescale of at 
most 10^ years set by sputtering, while the observed fraction of dust and AGN in early- 
type galaxies indicates that about 50% of early-type galaxies have dust at any time and 
this fraction approximately corresponds to the fraction with AGN. The main problem with 
a purely internal origin is that the dust creation from mass loss and the destruction from 
sputtering are continuous processes and all early- type galaxies should have some dust. The 
internal origin model therefore requires some more cataclysmic event to destroy the dust 
on approximately the same timescale required to reform the observed dust mass. In this 
case, the central AGN provides the most ready supply of the necessary energy to unbind 
or destroy the dust. The external origin, in contrast, requires continual dust replenishment 
via mergers with gas and dust-rich galaxies. Given the high fraction of carly-typc galaxies 
with dust, the frequency of these mergers must be comparable to the dust destruction time. 
In this case, the estimate of 10^ years appears to require an incredibly high merger rate 
with gas-rich dwarf galaxies. If this is the case, the dusty or active early-type galaxies may 
be preferentially in richer environments or have more companions. We checked this with 
measurements of the average local galaxy density and average projected angular separation 
to the nearest neighbor from (Ho et al. 1997a) and did not find evidence for higher local 
density or closer companions in the active or dusty samples. 



5.2. Nuclear Stelleir Disks in Ecirly-Type Galaxies 

Approximately a third of the inactive, early-type galaxies have structures that appear 
to be nuclear stellar disks. As noted in § 4, these are likely to be stellar disks, rather than 
nuclear stellar bars, because they have the same position angle as their host galaxy. The 
true incidence of stellar disks in our sample is likely to be higher than observed and may 
be 100% because most of these disks have very high axis ratios and we are most sensitive 
to stellar disks seen edge on. Many previous studies have found nuclear stellar disks in the 
central regions of elliptical and lenticular galaxies. Ravindranath et al. (2001) found such 
disks in 21% of their 33 early-type galaxies, Laucr et al. (2005) found that nearly all of the 
power-law galaxies in their sample have stellar disks, and Ferrarese et al. (2006) find stellar 
disks in 13% of a sample of 100 early-type galaxies in the Virgo Cluster. Nuclear stellar 
disks are therefore a common structure, although only a few have been studied in great 
detail (Pizzella et al. 2002; MoreUi et al. 2004; Kormendy et al. 2005). 

The origin of nuclear stellar disks has been a matter of recent debate, as these disks are 
not simple extensions of the large-scale disks to the centers of the galaxies. Krajnovic & Jaffe 
(2004) argue that the nuclear disks could be the result of the infall of mass to the center of 
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the galaxy driven by the secular evolution of a bar, galaxy mergers, or both. These authors 
also point out that none of the previously investigated nuclear stellar disk galaxies has an 
active nucleus, in agreement with our findings, although they do harbor a 10^~^ Mq SMBH, 
and argue that this makes them descendants of quasars that spent their fuel and turned off 
the central engine. Some known nuclear disks do exhibit strong evidence for an external 
origin, such as that in NGC 4698 (Pizzella et al. 2002), which is geometrically decoupled 
from the host galaxy. 

The ellipticity distribution of the five elliptical host galaxies with nuclear stellar disks 
may also have implications for the origin and dctcctability of these disks. As noted above, 
four of the five elliptical galaxies with nuclear stellar disks have very large axis ratios (three 
are classified E5 and one is classified 'E6?'). Although this is a small sample, such large 
axis ratios are rare and suggests that nuclear stellar disks could preferentially exist in the 
subset of intrinsically flattened ellipticals and not just observed in the subset of the nuclear 
disks with a favorable viewing angle. As lower- luminosity ellipticals are more likely to be 
oblate or flattened (Anthony & Rydcn 2005) and tend to be faster rotators (Cappellari et al. 
2006), the incidence of nuclear stellar disks should be reinvestigated as a function of elliptical 
galaxy properties, such as luminosity. 



5.3. Evolution and timescales 

The nuclear stellar disks observed in early- type, inactive galaxies may share a common 
origin with the dusty disks observed in early-type AGN hosts. Once a nuclear stellar disk is 
formed in a galaxy, it provides a natural preferred plane of rotation for new cold gas and dust. 
New material, either infalling or formed in situ would likely settle onto the disk in several 
rotation periods. While the dust settles, it would likely have a relatively chaotic appearance 
at first, although eventually it would form a relatively regular disk of dust and may host star 
formation that would increase the mass of the stellar disk. Kormendy et al. (2005) observe 
cospatial nuclear stellar and dust disks in the dwarf eUiptical NGC 4486A and suggest that 
these disks are toward the end stage of this evolutionary scenario, namely the observed 
stellar disk is forming from accreted gas and dust. This scenario is also supported by the 
recent study of Ferrarese et al. (2006). These authors find evidence of recent star-formation 
associated with dusty disks but not with the irregular, unsettled dust lanes. They suggest a 
similar evolutionary scenario for dust settling in which dust is acquired externally, form large 
scale dust lanes that gradually evolve to dusty disks with associated star-formation, and then 
finally collapse to nuclear dusty disks associated with stellar disks. As star formation ceases 
and their colors become similar to the host galaxy, the disks then become more difficult to 
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observe. 

Additional support for the evolutionary scenario above is the kinematic evidence for 
inflow along nuclear spirals in at least one case: NGC 1097. Its nuclear dusty spiral has been 
revealed by near-IR observations (Prieto et al. 2005) and by structure maps in Fathi et al. 
(2006). Using the Integral- Field Unit of the Gemini Multi- Object Spectrograph, the latter 
authors have found streaming motions along the nuclear spiral arms with inward velocities 
of up to 50kms~^ in the Ha emitting gas throughout the nuclear region. Another relevant 
finding on this galaxy is the young obscured starburst discovered by Storchi-Bergmann et 
al. (2005) very close to the nucleus (within ~10pc), in agreement with the suggestion that 
inflowing gas and dust gives birth to stars in the nuclear dusty spiral or disk. A key point in 
testing this evolutionary scenario is the evaluation of the life cycles of the dusty and stellar 
disks. In the case of NGC 1097, the velocity observed for the streaming motions along the 
nuclear spiral allows an estimate of a few Myr for the gas to flow from a few hundred parsecs 
to the nucleus. This is also consistent with the above estimates that the dust may survive 
for on order several dynamical timescales. 

If substantial star formation does add to stellar nuclear disks each activity cycle, then 
there should be some stars associated with the disks younger than the characteristic or 
episodic timcscale of the dust replenishment and AGN lifetime, or younger than 10^ years. 
While we do not have color information for our sample, a number of stellar disks arc indeed 
found to be bhic (Krajnovic & Jaffc 2004), although many authors also argue that the colors 
arc the same as those of the bulge (Erwin & Sparkc 2002). Relatively old ages - from 6 
to 15 Gyr - have been reported from studies using broad-band colors and spectral indices 
(Krajnovic & Jaffe 2004; MoreUi et al. 2004). Nevertheless, a close look at these studies 
reveals that the methods used to date the disk stellar population are not very sensitive 
to the presence of a small young component in the middle of a luminous old bulge. In 
particular, the spectral indices are measured in a restricted interval (5000 < A < 5400A) and 
this wavelength range is not very sensitive to a small population of young stars and a blue 
spectral index should be used instead. The same applies for the broad-band colors, which 
usually do not include images in the blue. Estimates of the total cold gas mass also suggest 
that the number of new stars formed from the observed dust is likely to be small compared 
to the mass of the stellar disks. Detailed stellar population studies of nuclear stellar disks, 
combined with measurements of the star formation rates in the centers of dusty early-type 
galaxies, could test the consistency of this simple evolutionary model. 
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6. Conclusions 

We have used archival HST WFPC2 images of 34 active and 34 inactive, early-type 
galaxies (T < 0) to investigate if activity is correlated with the presence of circumnuclear 
dust. We found a strong correlation between the presence of circumnuclear (hundreds of 
parsecs scales) dust and nuclear activity: all early-type AGN hosts have circumnuclear dust, 
while only 26% (nine of 34) early-type, inactive hosts have circumnuclear dust. This in- 
dicates that the presence of circumnuclear dust is a requirement, although not a sufficient 
condition, of black hole accretion in early-type galaxies. The circumnuclear dust in all of 
these early-type galaxies is typically confined to the central kiloparsec and is not observed 
on larger scales. The morphology of the circumnuclear dust is also typically complex, al- 
though a minority of galaxies (< 25%) show well-defined dusty disks with tightly wound 
spiral structure. Of the inactive galaxies, 38% (13) have structures that appear to be nu- 
clear stellar disks. As these structures appear to share the same position angle as the host 
galaxy, these structures are more likely to be stellar disks than nuclear stellar bars. We 
also investigated a sample of 31 active and 31 inactive late-type galaxies (T > 0) and found 
that all of the galaxies show evidence for circumnuclear dust, independent of the presence of 
nuclear accretion. 

While our data demonstrate a clear connection between dust and AGN in early-type 
galaxies, the origin of the dust and therefore the nature of the AGN fuehng mechanism 
remains unclear. We considered both internal and external origins for the circumnuclear 
dust, but both present significant challenges. Internal dust creation via stellar mass loss 
appears to be inconsistent with the absence of dust in approximately 50% of all early-type 
galaxies. An external origin, in contrast, appears to require frequent mergers of small gas 
and dust-rich galaxies. The timescale of these mergers is set by the settling and destruction 
times of the dust, which may be as short as a few dynamical times for these kiloparsec 
or smaller structures, or 10^ years. The large observed fraction of early-type galaxies with 
nuclear stellar disks, the even larger fraction of early-type galaxies with dusty disks, and the 
relatively short lifetime of the dust, appears to require continual growth and replenishment 
of nuclear stellar disks. Our observations demonstrate that this process is also intimately 
linked with some accretion onto the central, supermassive black hole. Observations of the 
dust kinematics and mass, along with improved models of the dust settling and destruction 
at the centers of early-type galaxies, could resolve the origin of the dust and nuclear stellar 
disks, and reveal the fueling mechanism for these low-luminosity AGN. 

We would like to thank Barbara Ryden for helpful discussions. RSL, TSB, MFS ac- 
knowledge support from the Brazilian institutions CNPq, CAPES and FAPERGS. Support 
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Table 1. HST data and general properties of the active sample 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

NGC WFPC2 Filter Exposure T D M^^ i Activity Pair Dust Disk Comments 

Number Camera (sec) (Mpc) (mag) (deg) Type Number (kpc) (kpc) 
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EARLY-TYPE EXTENDED SAMPLE 
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Table 1 — Continued 
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Note. — (3) Pass-band HST filter used in observation; (5—9) Galaxy T type, distance in Mpc, B-band absolute magnitude, inclination, 
and activity type from (Ho et al. 1997a), S indicates a Scyfert and L indicates a LINER; (10) NGC number of the control sample pair; 
(11) Main innermost dust structure semi-major axis size in kpc, a X sign means some dust structure is present but is difficult to be 
measured; (12) Stellar disk radial extent in kpc; (13) Comments: b = bar, d = dust disk, i = irregular dust, 1 = dust lanes or filaments, 
r = dust ring, s = dust spiral, w = weaJi feature, ? = uncertain classification. 
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Table 2. HST data and general properties of control sample 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 


WFPC2 


Filter 


Exposure 


T 


D 




i 


Activity 


Pair 


Dust 


Disk 


Comments 


Number 


Camera 




(sec) 




(Mpc) 


(mag) 


(deg) 


Type 


Number 


(kpc) 


(kpc) 





EARLY-TYPE MATCHED SAMPLE 



0821 


PC 


F555W 


1023 


PC 


F555W 


2300 


PC 


F555W 


2950 


PC 


F555W 


3377 


PC 


F555W 


3384 


PC 


F555W 


3640 


PC 


F555W 


4026 


PC 


F555W 


4291 


PC 


F555W 


4365 


PC 


F555W 


4371 


PC 


F606W 


4379 


PC 


F555W 


4382 


PC 


F555W 


4406 


PC 


F555W 


4469 


PC 


F606W 


4526 


PC 


F555W 


4570 


PC 


F555W 


4578 


PC 


F606W 


4612 


PC 


F606W 


4621 


PC 


F555W 


4649 


PC 


F555W 


4694 


PC 


F606W 


5308 


PC 


F555W 


5557 


PC 


F555W 


5576 


PC 


F555W 


7619 


PC 


F555W 


EARLY-TYPE EXTENDED SP 


0221 


PC 


F555W 


0507 


PC 


F555W 


3115 


PC 


F555W 


3412 


PC 


F606W 


3610 


PC 


F555W 


4245 


PC 


F606W 


4405 


PC 


F606W 


7457 


PC 


F555W 


LATE-TYPE 


MATCHED SAM 


0278 


WF 


F606W 


0864 


PC 


F606W 


0972 


PC 


F606W 


2146 


PC 


F606W 


2339 


PC 


F606W 


2748 


PC 


F606W 


2775 


PC 


F606W 



304 


-5 


23.2 


-20.11 




420 


-3 


10.5 


-20.03 


73 


304 


-2 


31 


-20.69 


44 


500 


-2 


23.3 


-20.04 


50 


160 


-5 


8.1 


-18.47 




380 


-3 


8.1 


-18.79 


65 


400 


-5 


24.2 


-20.73 




400 


-2 


17 


-19.56 


81 


304 


-5 


29.4 


-20.09 




1000 


-5 


16.8 


-20.64 




80 


-1 


16.8 


-19.51 


57 


160 


-2.5 


16.8 


-18.6 


32 


568 


-1 


16.8 


-21.14 


40 


500 


-5 


16.8 


-21.39 




160 





16.8 


-19.31 


73 


520 





16.8 


-19.31 


74 


200 


-2 


16.8 


-19.33 


76 


80 


-2 


16.8 


-18.96 


43 


80 


-2 


16.8 


-18.75 


38 


330 


-5 


16.8 


-20.6 




1050 


-5 


16.8 


-21.43 




500 


-2 


16.8 


-19.08 


63 


400 


-3 


32.4 


-20.13 




500 


-5 


42.6 


-21.17 




1600 


-5 


26.4 


-20.34 




1100 


-5 


50.7 


-21.6 





26 


-6 


0.7 


-15.51 




567 


-2 


65.7 


-21.96 




280 


-3 


6.7 


-19.39 


73 


80 


-2 


8.1 


-18.2 


57 


2400 


-5 


29.2 


-20.79 




160 





9.7 


-17.92 


41 


160 





31.5 


-19.63 


51 


1050 


-3 


12.3 


-18.69 


59 


320 


3 


11.8 


-19.69 


18 


560 


5 


20 


-20.25 


41 


600 


2 


21.4 


-20.17 


61 


560 


2 


17.2 


-20.6 


57 


600 


4 


30.9 


-20.97 


41 


600 


4 


23.8 


-20.29 


70 


300 


2 


17 


-20.34 


40 



A 


4636 




0.2 




A 


4111 




0.069 




A 


3998 








A 


3516 






b 


A 


4278 


0.26 


0.071 


1 


A 


2787 




0.034 




A 


4261 








A 


4036 


X 


0.071 


i,l 


A 


5077 








A 


4486 








A 


3414 


0.05 




d 


A 


5273 








A 


2655 








A 


2768 


0.02 




w 


H 


4293 


0.57 




i,l 


H 


3607 


1.21 






A 


4550 




0.119 




A 

A 


7743 








A 

A 


4143 








A 

A 


3945 




0.274 




A 

A 


4374 








H 


4138 


X 




i 


A 

A 


2911 




0.993 




A 


4589 








A 


3226 








A 


0315 




0.313 




A 


404 




0.004 




A 


5548 








A 


1358 




0.067 




A 


3166 




0.063 




A 


6340 




0.895 




H 


4203 


0.19 




r 


H 


3884 


X 






A 


4477 








H 


4394 


X 






H 


4639 


X 






H 


2639 


X 






H 


5377 


X 






H 


1667 


X 






H 


4501 


X 






A 


7217 


X 
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Table 2 — Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 


WFPC2 


Filter 


Exposure 


T 


D 




i 


Activity 


Pair 


Dust 


Disk 


Comments 


Number 


Camera 




(sec) 




(Mpc) 


(mag) 


(deg) 


Type 


Number 


(kpc) 


(kpc) 




2782 


PC 


F555W 


460 


1 


37.3 


-20.85 


43 


H 


2273 


X 






2903 


PC 


F606W 


560 


4 


6.3 


-19.89 


63 


H 


4258 


X 






3310 


PC 


F606W 


500 


4 


18.7 


-20.41 


40 


H 


1068 


X 






3344 


PC 


F606W 


160 


4 


6.1 


-18.43 


24 


H 


1058 


X 






3351 


PC 


F606W 


560 


3 


8.1 


-19.28 


48 


H 


3368 


X 






3423 


PC 


F606W 


160 


6 


10.9 


-18.81 


32 


H 


3486 


X 






3504 


PC 


F606W 


500 


2 


26.5 


-20.61 


40 


H 


6500 


X 






3556 


PC 


F606W 


160 


6 


14.1 


-20.92 


80 


H 


3079 


X 






4041 


WF 


F606W 


320 


4 


22.7 


-20.02 


21 


H 


3507 


X 






4102 


PC 


F606W 


600 


3 


17 


-19.54 


56 


H 


4450 


X 






4254 


PC 


F606W 


560 


5 


16.8 


-21.03 


30 


H 


6951 


X 






4274 


PC 


F555W 


296 


2 


9.7 


-19.35 


71 


H 


3718 


X 






4380 


PC 


F606W 


560 


3 


16.8 


-19.06 


58 


H 


4698 


X 






4448 


PC 


F606W 


160 


2 


9.7 


-18.65 


71 


H 


4772 


X 






4559 


PC 


F606W 


160 


6 


9.7 


-20.17 


68 


H 


2033 


X 






4567 


PC 


F606W 


160 


4 


16.8 


-19.34 


48 


H 


4651 


X 






4654 


PC 


F606W 


160 


6 


16.8 


-20.38 


56 


H 


7479 


X 






4800 


PC 


F606W 


160 


3 


15.2 


-18.78 


43 


H 


3982 


X 






4845 


PC 


F606W 


160 


2 


15.6 


-19.55 


79 


H 


4388 


X 






5248 


PC 


F517M 


900 


4 


22.7 


-21.15 


44 


H 


5985 


X 






5383 


PC 


F606W 


560 


3 


37.8 


-20.94 


32 


H 


3642 


X 






5806 


PC 


F606W 


600 


3 


28.5 


-20.43 


61 


H 


5566 


X 






6217 


PC 


F606W 


500 


4 


23.9 


-20.23 


34 


H 


4579 


X 






6412 


PC 


F606W 


560 


5 


23.5 


-19.78 


30 


H 


4051 


X 







Note. — (3) Pass-band HST filter used in observation; (5—9) Galaxy T type, distance in Mpc, B-band absolute magnitude, 
inclination, and activity type from (Ho ct al. 1997a), S indicates a Scyfcrt and L indicates a LINER; (10) NCC number of the control 
sample pair; (11) Main innermost dust structure semi-major axis size in kpc, a X sign means some dust structure is present but is 
difficult to be measured; (12) Stellar disk radial extent in kpc; (13) Comments: b = bar, d = dust disk, i = irregular dust, 1 = dust 
lanes or filaments, r = dust ring, s = dust spiral, w = weak feature, ? = uncertain classification. 
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Fig. 1. — Distributions of plate scales, Hubble types, absolute B magnitudes, and inclinations 
for all active and control galaxies of the matched sample. 
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Fig. 2. — Distributions of plate scales, Hubble types, absolute B magnitudes, and inclinations 
for all active and control galaxies of the extended sample. 
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Fig. 3. — Structure maps for the matched sample of early-type galaxies. Each image covers 
5% of the galaxy D25, and is rotated so that North is up and East is to the left. In each panel 
the left column displays the active galaxy and the right column the pair-matched control 
galaxy. 




Fig. 3. — continued 




Fig. 3.- 



continued 




Fig. 4. — Structure maps for the extended sample of early-type galaxies. Images are pre- 
sented as in Fig. 3. 




Fig. 5. — Structure maps for the matched sample of late- type galaxies. Images are presented 
as in Fig. 3. 




Fig. 5. — continued 




Fig. 5. — continued 
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Fig. 6. — Presence of dust structure in the matched (upper panel) and extended (lower 
panel) sample as a function of Hubble type for the active and control galaxies. 
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Fig. 7. — Presence of stellar disks in the matched (upper panel) and extended (lower panel) 
sample as a function of Hubble type for the active and control galaxies. 
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Fig. 8. — Stellar disk's radial extent in kpc for the matched (upper panel) and extended 
(lower panel) sample. 



